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Q ; ABSTRACT 

, We present new observations of the transiting system TrES-3 obtained from 2009 to 

c/2 ■ 2011 at several observatories. The orbital parameters of the system were re- determined 

and a new linear ephemeris was calculated. The best quality light curve was used for 
light curve analysis, and other datasets were used to determine mid-transit times (Tc) 
and study transit time variation (TTV). For planet parameter determination we used 
\ two independent codes and finally, we concluded that our parameters are in agreement 

■ with previous studies. Based on our observations, we determined 14 mid-transit times. 

C*~) ' Together with published Tc we found that the timing residuals showed no significant 

IO , deviation from the linear ephemeris. We concluded that a periodic TTV signal with 

an amplitude greater than 1 minute over a 4- year time span seems to be unlikely. Our 
(f-) I analysis of an upper mass limit allows us to exclude an additional Earth-mass planet 

c^) , close to inner 3:1, 2:1, and 5:3 and outer 3:5, 1:2, and 1:3 mean-motion resonances. 

(^T) ■ Using the long-term integration and applying the method of maximum eccentricity, the 

region from 0.015 au to 0.05 au was found unstable and the region beyond the 0.05 au 
was found to have a chaotic behaviour and its depletion increases with increasing 
values of the initial eccentricity as well as inclination. 
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1 INTRODUCTION 

Since the f irst discovery of a tran siting planet around 
HD 209458 (|Charbonneau et al.llioool 'l. 235 transiting extra- 
solar systems have already been confirmed up to December 
4th, 20120 Whilst transiting exoplanets offer unique scien- 
tific possibilities, their study involves several complications. 
In general, it is impossible to measure the mass and radius 
of a planet based on a dataset obtained with one observa- 



* Partly based on observations made at the Centro Astronomico 
Hispano Aleman (CAHA), operated jointly by the Max-Planck 
Institut fur Astronomie and the Instituto de Astroffsica de An- 
daluci'a (CSIC). 
f vanko@ta3.sk 
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tional technique. Transit light curves allow us to determine 
just the relative size of a star and planet, the orbital incli- 
nation and the stellar limb-darkening coefficients. By com- 
bining this with radial- velocity measurements, the observa- 
tions offer the opportunity to measure the precise stellar and 
planetary parameters. In order to obtain such parameters, 
some constraints are needed, and are usually provided by 
forcing the prope rties of the host s tars to match theoreti- 
cal expectations l|Southworthl l201(ih . Significant uncertain- 
ties remain in the stellar mass and radius determinations 
of many systems. In some cases, this is due to poorly de- 
termined photospheric properties (i.e. effective temperature 
and metallicity), and in o ther cases due to a la ck of an accu- 
rate luminosity estimate (|Sozzetti et al.ll2009l h In addition, 
the different methods used for these determinations as well 
as different approaches toward systematic errors are leading 
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to rather inhomogeneous set of planet properties. Because of 
such inhomogenities, recently a few papers were published 
where authors re-analyzed a large subset of known transiting 
plan ets, and applied a uniform methodology to all systems 
(e.g. iTorres et al.ll200Sl ; ISouthworthll2010l . l2012h . 

In this paper we focus on the transiting system TrES-3. 
The system consists of a nearby G-type dwarf and a massive 
hot Jup iter with an orbita l perio d of 1.3 days. It was discov- 
ered by lO 'Donovan et ah I J2007I) and also detected by the 
SuperWASP survey (|Collier Cameron et"al] 120071 ). Later, 
ISozzetti et alj (|2009l ) presented new spectroscopic and pho- 
tometric observations of the host star. A detailed abundance 
analysis based on high- resolution spectra yields [Fe/H] = 
-0.19 ± 0.08, r e a = 5650 ± 75 K and log g = 4.0 ± 0.1. The 
spectroscopic orbital solution was im proved with new radia l 
velocity measurements obtained by ISozzetti et al.l (|2009h . 
Moreover, these authors redetermined the stellar parame- 
ters (i.e. M* = 0.928±£ 



| Sozzetti et all 120091: IGibson et al l l2009l ; IColon et al.ll2010l; 
Southworth 20 10]. 12011c iLee et al.ll201ll; ICh ristianse n et al 
201ll;ISada et al.ll2012l; (Turner et al.ll2013l; IKundurthv et al 



20131 ). (ii) based on the obtained transits, we will determine 



M and R* = 0.829Z 



Re 



and finally, the new values of the planetary mass and ra- 
dius were determined (see Tab [3J. They also studied the 
transit timing variations (TTVs) of TrES-3 and noted sig- 
nificant outli e rs fro m a constant period. In the same year, 
IGibson et all 12009) presented the follow-up transit photom- 
etry. It consisted of nine transits of TrES-3, taken as part 
of transit timing program using the RISE instrument on 
the Liverpool Telescope. These t ransits, together wit h eight 
transit times published before (|Sozzetti et al.l [2009T ) . were 
used to place upper mass limit as a function of the period 
ratio of a potential perturbing planet and transiting planet. 
It was shown that timing residuals are sufficiently sensitive 
to probe sub-Earth mass planet in both interior and ex- 
terior 2:1 resonances, assuming that the additional planet 
is in an initially circular orbit. IChristiansen et ail (|201ll ) 
has observed TrES-3 as a part of the NASA EPOXI Mis- 
sion of Opportunity. They detected a long-term variability 
in the TrES-3 light curve, which may be due to star spots. 
They also confirmed that the planet ary atmosphere does not 
have a temperature inversion. Later. iTurner et al.l (|2013T ) ob- 
served nine primary transits of the hot Jupiter TrES-3b in 
several optical and near-UV photometric bands from June 
2009 to April 2012 in an attempt to detect its magnetic 
field. Authors determined an upper limit of TrES-3b's mag- 
netic field strength between 0.013 and 1.3 G using a tim- 
ing difference of 138 s derived from the Nyquist-Shannon 
sampling theorem. They also presented a refinement of the 
physical parameters of TrES-3b, an updated ephemeris and 
its first published near-UV light curve. The near-UV plan- 
etary radius of R p = 1.386±^g Rj was also determined. 
This v alue is consistent with th e planet's optical radius. Re- 
cently, IKundurthv et all (|2013l ) observed eleven transits of 
TrES-3b over a two year period in order to constrain system 
parameters and look for transit timing and depth variations. 
They also estimated the system parameters for TrES-3b and 
found consistency with previous estimates. Their analysis of 
the transit timing data show no evidence for transit tim- 
ing variations and timing measurements are able to rule out 
Super-Earth and Gas Giant companions in low order mean 
motion resonance with TrES-3b. 

The main aims of this study can be summarized in the 
following items: (i) determination of the system parameters 
for TrES-3b (two independent code s will be used) and com- 
parison with previous studies (i.e. lO'Donovan et all 120071 : 



the mid-transit times (Tc) and with following analysis of 
transit time variation (TTV) we will discuss possible pres- 
ence of a hypothetical additional planet (perturber). We will 
try to estimate its upper-mass limit as a function of orbital 
periods ratio of transiting planet and the hypothetical per- 
turber. (iii) Finally, using the long-term integration and ap- 
plying the method of maximum eccentricity we will search 
for stability of regions inside the TrES-3b planet in context 
of additional planet (s). 

The remainder of this paper is organized as follows. In 
the Section 2, we describe observations and data reduction 
pipelines used to produce the light curves. Section 3 presents 
the methods for analysis of transit light curves as well as dis- 
cussion and comparison of the parameters of TrES-3 system. 
Section 4 and 5 are devoted to TTV and long-term stability 
of the system, respectively. Finally, in Section 6 we summa- 
rize and discuss our results. 



2 OBSERVATIONS AND DATA REDUCTION 

We obtained our data using several telescopes with different 
instruments. This allowed us to obtain many light curves 
since this strategy can effectively cope with the weather 
problems. On the other hand, this approach results in rather 
heterogeneous data. We used most of the data in average 
quality for the TTV analysis which is not very demanding 
on homogeneity of the data. Only the best quality light curve 
was used for the planet parameter determination. 

Observations used in this paper were carried out at 
the several observatories in Slovakia (Stara Lesna Obser- 
vatory; 49°09' 10"N, 20°17' 28" E), Poland (Piwnice Obser- 
vatory; 53°05' 43"N, 18° 13' 46" E), Germany (Grossschwab- 
hausen Observatory; 50°55' 44"N, 11°29' 03"E; Volksstern- 
warte Kirchheim Observatory, 50°55'44"N, 11°29' 03"E and 
Michael Adrian Observatory, 49°55'27"N, 08°24' 33"E) and 
Spain (Calar Alto Observatory; 37°13'25"N, 02°32' 46" E). 
We collected 14 transit light curves obtained between May 
2009 and September 2011. The transits on May 12, 2009 
and August 20, 2010 were observed simultaneously at two 
different observatories. The telescope diameters of 0.5 to 2.2 
m allowed us to obtain photometry with 1.2 - 7.8 mmag 
precision, depending on observing conditions. Observations 
generally started ~ 1 hour before the expected beginning of 
a transit and ended ~ 1 hour after the event. Unfortunately, 
weather conditions and schedule constraints meant that we 
were not able to fit this scheme in all cases. 

All instruments are equipped with CCD cameras with 
the Johnson-Cousins (UBVRcIc) standard filter system. 
The information from individual observatories and instru- 
ments as well as the summary of observing runs are given 
in Table [T] and Tabled The standard correction procedure 
(bias, dark and flat field correction) and subsequently aper- 
ture photometry was performed by IRAljf] and task chphot 



2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
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Table 1. Overview of the telescopes and instruments/detectors 
used to obtain photometry of TrES-3. FoV is the field of view of 
the instrument and Ntr is the number of observed transits. Ab- 
breviations of the observatories: Gl - Stara Lesna Observatory, 
GSH - Grofischwabhausen observing station of the Jena Uni- 
versity (CTK — Cassegrain Telesko p Kamera; STK — Schmidt 
Teleskop Kamera, see [Mugrauer] 2009; Mugraucr & B ertholdl 
2010), MA - Michael Adrian Observatory in Trebur, VK - 
Volkssternwarte Kirchheim Observatory (RCT - Ritchie Chretien 
Telescope), P - Piwnice Observatory and CA - Calar Alto Ob- 
servatory (RCF - Ritchie Chretien Focus). 



Obs. 


Telescope 


Detector 








CCD size 


FoV [arcmin] 


Gl 


Newton 


SBIG ST10-MXE 


5 






508/2500 


2184 x 1472, 6.8 fim 


20.4 X 


13.8 


MA 


Cassegrain 


SBIG STL-6303E 


2 






1200/9600 


3072 x 2048, 9 fim 


10 X 


7 


GSH 


CTK 


SITe TK1024 


1 






250/2250 


1024 x 1024, 24 fim 


37.7 x 


37.7 




STK 


E2V CCD42-10 


1 






600/1758 


2048 x 2048, 13.5 fim 


52.8 x 


52.8 


VK 


RCT 


STL-6303E 


2 






600/1800 


3072 x 2048, 9 fim 


71 x 


52 


P 


Cassegrain 


SBIG STL-1001 


1 






600/13500 


1024 X 1024, 24/xm 


11.8 x 


11.8 


CA 


RCF 


SITe CCD 


2 






2200/17037 


2048 X 2048, 24/^m 


18.1 x 


18.1 



Table 2. Summary of the observing runs: Obs. - Observatory 
according to Table 1, N exp - number of useful exposures, t exp - 
exposure times. The dates are given for the beginning of nights. 



Obs. 


Date 


Filter 


exp 


texp (s 


Gl 


2009 May 12 


R 


319 


40 




2009 Aug 01 


R 


345 


45 




2010 Apr 27 


R 


180 


35 




2010 Jun 30 


R 


238 


40 




2010 Aug 07 


R 


168 


35 


MA 


2010 July 13 


R 


349 


25 




2010 Aug 20 


R 


296 


20 


GSH (CTK) 


2009 May 25 


1 


138 


80 


(STK) 


2011 Mar 22 


R 


131 


60 


VK 


2009 Aug 14 


Clear 


103 


120 




2010 Aug 20 


Luminance 


320 


30 


P 


2009 May 12 


R 


120 


55 


CA 


2010 Sep 06 


R 


158 


30-35 




2011 Sep 12 


R 


128 


45 



For the calculation of synthetic light curves we used 
two independe nt approaches: the fir st one is based on the 
routines from (jMandel fc Agoll I200I. ) and the Monte Carlo 
simulations (de scribed in Section 13.1 [>. the second one uses 
JKTEB0P code (|Southworth. Maxted fc Smallevl 12004 ) (see 
Section I3.2[) . 



l|Raetz et alj 120091) (GSH and VK), C-munipack packag«0 
(Gl) and Mira_Pro_70 (MA). Data from remaining tele- 
scopes (P and CA) were reduced with the software pipeline 
develop ed for the Semi-Automatic Variability Search sky 
survey (jNiedzielski. Macieiewski fc Czard 120051 ) . To gener- 
ate an artificial comparison star, at least 20-30 per cent of 
stars with the lowest light-curve scatter were selected itera- 
tively from the field s tars brighter than 2 .5-3 mag below th e 
saturation level (e.g. iBroeg et al.ll2005l : iRaetz et aHl200sh . 
To measure instrumental magnitudes, various aperture radii 
were used. The aperture which was found to produce light 
curve with the smallest overall scatter was applied to gen- 
erate final light curve. The linear trend in the out-of-transit 
parts was also removed. 



3 LIGHT CURVE ANALYSIS 

Our photometric observations of the TrES-3 system consist 
of data from different instruments and are of different photo- 
metric quality. For the purpose of radius determination, we 
decided to analyse only a light curve with the lowest scatter 
and the best rms = 1.2 mmag. We choose data obtained at 
Calar Alto observatory on September 06, 2010 (see Figs. 1 
and 2). We first refined the light-curve system parameters 
and subsequently determined the individual times of transit 
as described below in order to improve ephemeris. 

Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

3 http: / / c-munipack.sourceforge.net/ 

4 http:/ /www. mirametrics.com/mira_pro. htm 



3.1 SOLUTION 1 

First we used the downhill simplex mi nimization proce - 
dure (implemented in routine AMOEBA. IPress et aUll992l ) 
to determine 4 system parameters R p /R, (planet to star 
radius ratio), i (inclination), Tc (mid-transit time) and 
Rt/a (star radius to semi-major axis ratio). The model 
light curve itself was compu ted via the analytic expressions 
from lMandel fc Agoll (|2002l ). The quadratic limb darkening 
law was assumed and corresponding limb da rkening coeffi- 
cients a, C2 were linearly interpolat ed from IClaretl (|2000|) 
assuming the stellar parameters from lSozzetti et alj l|2009() : 
T cff = 5 650 K, log (g) = 4.4 and [Fe/H] = -0.19. As a 
goodness of the fit estimator we used the \ 2 function: 

JV , 

where rrii is the model value and di is the measured value of 
the flux, Oi is the uncertainty of the i th measurement and 
the sum is taken over all measurements. The orbital period 
and the limb darkening coefficients were fixed through the 
minimization procedure. The transit duration To was deter- 
mined assuming the semi-major axis a = 0.022821q'qoo4o au 
l|Sozzetti et al.1 [20091). The parameters and the correspond- 
ing light curve for which we found the minimum value of x 2 
function, are in Table 3 and Figure 2, respectively (SOLU- 
TION 1). 

To estimate the uncertainties of the calculated tran- 
sit para meters, we employ ed the Monte Carlo simulation 
method l|Press et alJll992h . We produced 10 000 synthetic 
data sets with the same probability distribution as the resid- 
uals of the fit in Figure [2] From each synthetic data set 
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Figure 1. The light curves of the system TrES-3 obtained at individual observatories between May 2009 and September 2011. The best 
light curve (obtained at Calar Alto, Sep 2010) was used to determine system parameters. The data from other observatories were used 
for Tq determination and TTV analysis. The typical error bar is plotted next to each light curve. 
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obtained by this way we estimated the synthetic transit pa- 
rameters. The minimum Xmin value corresponding to each 
set of Monte Carlo parameters was calculated as: 



2 

Xmin 



N 

E 

i=l 



(2) 



where m, is the original best fit model value and Si is the 
'Monte-Carlo simulated' value. Figure [3] shows the depen- 
dence of the parameters Rp/R* and i on the reduced x?- 
This quantity is defined as: 



2 Xmin 

Xr 



N - M' 



(3) 



where N is the number of data points and M is the number 
of fitted parameters. To fully understand the errors of the 
system parameters we constructed confidence intervals in 2- 
dimensional space. Figure 2] depicts the confidence regions 
for 2 parameters (R p /R, vs. i) as a projection of the original 
4-dimensional region. The gray-colored data points stand 
for the la, 2a and 3a region with corresponding value of 
Ax 2 = Xmin - Xm = 4.72, 9.7 and 16.3, respectively (xL is 
the x 2 of the original best fit model value) . From the shape 
of the dependence it could be seen that these two parameters 
correlate (see Section T3.3p . 

Finally, we took into account the uncertainty of the 
stellar mass and semi-major axis according to the simple 
error propagation rule. The results from the first analysis 
are shown in the Table [3] as SOLUTION 1. 
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3.2 SOLUTION 2 (JKTEBOP code) 

We h ave modelled the light curve using the JKTEBOrQ 
code (|Southworthll201ll ) as well. JKTEBOP grew out of the 
original EBOP program written for eclipsing binary systems 
IIEtzej | l98ll; IPopper fe Etzej |l98l|) via implementing the 
NDE l|Nelson fc Davisll 19721 ) model. JKTEBOP uses biaxial 
spheroids to model the component stars (or star and planet) 
and performs a numerical integration in concentric annuli 
over the surface of each body to obtain the flux coming 
from the system. This feature of the code allows us to avoid 
small and spherical planet approximations which are used 
in ana lytic light-curve generators based on lMandel fc Agoll 
(2002), and hence to derive planets oblateness. A model is 
fitted to the data by the Levenberg-Marquardt least-square 
procedure. The code converges rapidly toward a reliable so- 
lution and diminishes th e correlation between fitted param- 
eters jSouthworthll20Qg| ). 

The main parameters of a JKTEBOP fit are the orbital 
inclination i, and the fractional radii of the host star and 
planet, ta and rv The fractional radii are defined as: 



a 



R 



Tb = 



(4) 



where R, and R p are the stellar and planetary radii and 
a is the orbital semi-major axis. Parameters va and rb 
correspond to radii of spheres of the same volume as the 



5 JKTEBOP is written in FORTRAN77 

and the source code is available at 
http://www.astro.keele.ac.uk/jkt/codes/jktebop.html 



Figure 2. Top: The light curve obtained at Calar Alto on Septem- 
ber 06, 2010, and the best fit corresponding by SOLUTION 1. 
Middle: Residuals from the best fit mentioned above. Bottom: 
Difference from the fits corresponding by SOLUTION 1 and SO- 
LUTION 2. 



biaxial spheroids. In JKTEBOP the fractional radii are 
reparametrized as their sum and ratio: 



f A + r b , 



n 

TA 



Rp 



(5) 



because these are only weakly correlated with each other 
l|Southworthll2008h . The directly fitted orbital inclination, i, 
allows the transit impact parameter b = cos i to be cal- 
culated. T he initial value s of pa rameters listed above were 
taken from ISozzetti et al.l l|2009h . Because of different qual- 
ity of data, the synthetic light curve was calculated only 
for the best light curve obtained at Calar Alto on Septem- 
ber 06, 2010 (the same l ike in SOLUTION 1). A value of 
a = 0.02282±g;°°°4o au (|Sozzetti et all 120091 ) was used in 
subsequent calculations. The best-fit model was used as a 
template and fitted to other light curves for which only 
the mid-transit time was allowed to vary. The resulting 
mid-t ransit times togeth e r with Tg obtaine d from litera- 
ture (jSozzetti et al. 20091: ICibson et al.|[2009l ) are analysed 
in details in the Section 4. The determined parameters ob- 
tained by JKTEBOP code are presented in Table [3] as SO- 
LUTION 2. 

The errors of derived parameters were determined in 
two ways for each combination of data set and adopted LD 
law (we have used the same coefficients like in the case of 
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Figure 3. Top: Reduced Xi representing global minimum solu- 
tion for calculated parameter R p /R*. Bottom: Reduced \? rep- 
resenting global minimum solution for calculated parameter i. 

SOLUTION 1). Firstly, we ran 1000 Monte Carlo (MC) sim- 
ulations, a spread range of a given parameter within 68.3% 
was taken as its error estimate. Secondly, the pr ayer-bead 
method (e.g. lDesert et al.ll20lil : IWinn et al.ll2009l ) was used 
to check whether red noise was present in our data. MC er- 
rors were found to be 2 - 3 times smaller than the values 
returned by the prayer bead method. This indicates that the 
light curve is affected not only by Poisson noise but also by 
additional correlated noise. Therefore, our prayer-bead error 
estimates were taken as our final errors. 



3.3 Light curve analysis results 

The resulting values of the parameters together with their 
uncertainties are given in Table [3] For comparison, the pa- 
rameters from previous studies are added there as well. Fig- 
ure[2]shows the resulting best fit obtained by SOLUTION 1. 
In order to compare both solutions, we also plotted differ- 
ences from the fit obtained by routines resulting from SO- 
LUTION 1 and SOLUTION 2. 

The final parameters are in good agreement with al- 
ready published values. The output from SOLUTION 1, in 
p articular R v and J?*, c orres pond a bit better to param eters 
of ISozzetti et all (|2009l ') and Ichristiansen et alj l|201lf) . The 
radii determined by JKTEBOP are also within the range of 
errors. Based on these two parameters we also inferred the 
critical value of inclination for total transit of TrES-3b as: 



Figure 4. Confidence region depicted as a projection of the 4- 
dimensional region into 2-dimensional parameter space (inclina- 
tion and ratio of the planet to star radius). Regions of la, 2<r and 
3<T corresponding to Ax = 4.72, 9.7 and 16.3, respectively, are 
marked with different gray colors. 

cosi = — - ^ =>i = 81.9°. (6) 
a a 

As listed in Table [3] all values of inclination including uncer- 
tainties are in agreement with critical inclination calculated 
above. This is an evidence of grazing transit of TrES-3b. 

We used MC simulations to produce of 10000 synthetic 
data sets with the same probability distribution as the resid- 
uals of the fit. From each synthetic data set we estimated the 
synthetic transit parameters. Using the results of the simu- 
lation, the dependence of the inclination versus the planet to 
star radius ratio is plotted in Figure [3] The figure demon- 
strates that there is a correlation between the parameters 
and that solution for the system TrES-3 is not unique and 
can be located in a relatively wide range (degeneracy of the 
parameters). This fact is also caused by grazing transit and 
by subsequent sensitivity of the solution to the determina- 
tion of LD coefficients. 

For determination of our mid-transit times, the best- 
fit model was used as a template and fitted to other light 
curves for which only the mid-transit time was allowed to 
vary. In order to have homogeneous analysis of TTV, we 
collected the data points of transit light curves published 
in previous studies (see Figure Q . Trans i t light curve data 
of ISozzetti et all l|2009l ) and IColon et al.1 l|2010h were avail- 
able on-line and the remaining tabulated data were obtained 
from mentioned authors by e-mail. We re-analysed all transit 
light curves and derived mid-transit times by the same pro- 
cedure (the code JKTEBOP) to get a homogeneous dataset 
for TTV analysis (see Section 0}. 

We also re-analysed all light curves under consideration 
with k — rb/rA as a free parameter and tried to search 
for any variation. All determinations of k were consistent 
within error bars with a mean value, so we did not detect 
any transit depth variation. We also saw no significant signal 
in periodograms. 

4 TRANSIT TIMING VARIATION 

Our new 14 mid-transit times and 42 redetermined lit- 
erature values were used to refine the transit ephemeris. 
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Table 3. Parameters of the extrasolar system TrES-3 from this work (This work: SOLUTION 1 and SOLUTION 2) compared with the 
results from the previous studies. R,/a is the star radius to semi-major axis ratio, Rp/R* is the planet to star radius ratio, i i s the 
inclination of the orbit, Trj is the transit duration assuming the semi-major axis of a = 0.02282^q'qqq 2 q au jSozzetti et al,lf2009l ) and 
P or b is the orbital period. The orbital period in this work was fixed during analysis. The errors of the orbital periods are in parenthesis. 



Source 


R*/a 


Rp/R* 


i 


T D 


Porb 








[°] 


[min] 


[days] 



O 'Donovan et al. f 20071 
Sozzetti et al. (2009) 
Gibson et al. (2009) 
Colon et al. (2010) 
Southworth (2010) 
Lee et al. (2011) 
Christiansen et al. (2011) 
Southworth (2011) 
Sada et al. (2012) 
Kundurthv et al. (2013) 
Solution.l 

Kundurthv et al. (2013) 
Solution_2 
Turner et al. 
This work 
Solution_l 
This work 
Solution_2 



0.1650 ± 0.0027 0.1660 ± 0.0024 
0.1655 ± 0.0020 



1687+ ' 0140 

U -- LOO '_0.0410 



1664+ 0011 
u.roo<i_ 001g 

1662+ 0046 

U.1DDZ_ 004g 



0.1603 ± 0.0042 
0.1664 ± 0.0204 0.1661 ± 0.0343 



82.15 ± 0.21 
81.85 ± 0.16 
81. 73 4 



,+0.13 
-0.04 



82.07 ± 0.17 
81.77 ± 0.14 
81.99 ± 0.30 
81.93 ± 0.13 



0.1675 ± 0.0008 0.1652 ± 0.0009 81.95 ± 0.06 



0.1698 ± 0.0014 0.1649 ± 0.0015 



0.1721 



+0.0054 
0.0052 



0.1693 



+0.0087 
0.0069 



0.1682 ± 0.0032 0.1644 ± 0.0047 



0.1696 



+0.0024 



0.1669 



+0.0027 



81.51 ± 0.14 
81.35lg;g 

81.86 ± 0.28 



79.92] 
83.77" 



O.liO 
1.15 
2.79 



81.9 ± 1.1 
77.9 ± 1.9 

81.30 ± 0.23 
79.20 ± 1.38 
79.08 ± 0.72 



1.30619(1) 
1.30618581(51) 
1.3061864(5) 

1.3061864(5) 
1.30618700(15) 
1.30618608(38) 
1.30618700(72) 

1.3061865(2) 

1.3062132(2) 

1.3062128(2) 
1.306 1854(1) 

1.306186 

1.306186 



BJD TDB -2450000 




U 

6 



Figure 5. Observation minus calculation O — C diagram for transit timing of TrES-3 b, plotted according to the new linear ephemeris. 
Open circles denote re-analysed mid-transit times from the literature. Filled symbols mark new mid-transit times reported in this paper. 



The mid-transit times were transformed from JD or HJD 
(based on UTC) into BJD (based on Barycentric Dy- 
namical Time - TDB) using the on-line converted by 
lEastman. Siverd fc GaudU |20ld ). As a result of fitting a lin- 
ear function of the epoch, we obtained the mid-transit time 
for the initial epoch T = 2454538.58144 ± 0.00007 BJDtdb 
and the orbital period P b = 1.30618599 ±0.00000023 d. The 
individual mid-transit errors were taken as weights. The lin- 
ear fit yields re duced y 2 of 2. 5 that is similar to a value of 
2.3 reported bv lGibson et al. | l|2009l ). These values, notice- 
able greater than 1 might suggest the existence of an addi- 
tional planet which perturbs the orbital motion of TrES-3 b 
jSozzetti et al.ll2009l ; iGibson et al1l2009l ) . 



http: / /astroutils. astronomy.ohio-state.edu/time/utc2bjd. html 



Results for new mid-transit times are shown in Ta- 
ble [4] The observed minus calculated (O-C) diagram, plot- 
ted in Fig. [SI shows no significant deviation from the linear 
ephemeris. All our data points deviated by less than 3a. We 
also searched for a periodicity that could be a sign of an ad- 
ditional body in the system. We genera ted a Lomb-Scargle 
periodogram (|Lomblll976l ; IScargld[l982T ) for the residuals in 
a frequency domain limited by the Nyquist frequency and 
found the highest peak at ~30 d and peak-to-peak ampli- 
tude of 70 ± 20 s. This period could coincide with a stellar 
rotation period which is roughly estimated to be ~28 d. Ex- 
amples of such TTV signals induced b y the stellar activit y 
are observed in the Kepler data (e.g. iMazeh et al.l 120131 ). 
However, the false alarm probability (calculated empirically 
by a bootstrap resampling method with 10 J trials) of the 
putative signal is disqualifying with a value of 18.2%. In ad- 
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Table 4. Results of transit timing. Obs. codes an observatory and instrument according to Table 1. To denotes the mid-transit times 
given as BJD (based on Barycentric Dynamical Time, TDB). Errors of mid-transit times are in days. The O — C values were calculated 
according to the new ephemeris. 



Date 




Obs. 


Epoch 


T (BJDtdb) 


To error 


O- 


-C(d) 


2009 


May 12 


P 


326 


2454964.39885 


0.00084 


+0 


00077 


2009 


May 12 


Gl 


326 


2454904. 39859 


0.00081 


+0 


00051 


2009 


May 25 


GSH (CTK) 


336 


2454977.4593 


0.0016 


-0 


0007 


2009 


Aug 01 


Gl 


388 


2455045.38091 


0.00045 


-0 


00070 


2009 


Aug 14 


VK 


398 


2455058.44418 


0.00077 


+0 


00071 


2010 


Apr 27 


Gl 


594 


2455314.45510 


0.00062 


-0 


00082 


2010 


Jun 30 


Gl 


643 


2455378.45937 


0.00031 


+0 


00034 


2010 


July 13 


MA 


653 


2455391.52067 


0.00017 


-0 


00022 


2010 


Aug 07 


Gl 


672 


2455416.3397 


0.0011 


+0 


0012 


2010 


Aug 20 


MA 


682 


2455429.40118 


0.00033 


+0 


00090 


2010 


Aug 20 


VK 


682 


2455429.3997 


0.0004 


-0 


0006 


2010 


Sep 06 


CA 


695 


2455446.38075 


0.00017 


+0 


00005 


2011 


Mar 22 


GSH (STK) 


846 


2455643.6145 


0.0005 


-0 


0003 


2011 


Sep 12 


CA 


979 


2455817.3372 


0.0014 


-0 


0003 



dition, the amplitude is close to the mean 1-a timing error of 
our observations. These findings allow us to conclude that a 
strictly periodic TTV signal with the amplitude greater than 
~1 minute ove r a 4-year tim e span seems to be unlikely. 

Following Gib son et alj (|2009T ). we put upper con- 
straints on a mass of a potential perturbing planet in the 
system with refined assumptions. We adopted a value of 1 
min for the maximal amplitude of the possible TTV signal. 
We also increased the sampling resolution to probe reso- 
nances other than i nner and out e r 1:2 commensurabilities 
studied in detail bv lGibson et~aH (|2009h ■ We simplified the 
three-body problem by assuming that the planetary system 
is coplanar and initial orbits of both planets are circular. The 
masses of both the star and the trans iting planet, as well a s 
its semi-major axis were taken from ISozzetti et al.l l|2009[) . 
The orbital period of the hypothetical perturber was varied 
in a range between 0.2 and 5 orbital periods of TrES-3 b. 
We produced 2800 synthetic O-C diagra ms based on ca l- 
culations done with the Mercury code (|Chamberslll999f ). 
We applied the Bulirsch-Stoer algorithm to integrate the 
equations of motion. 

The most important feature of the Bulirsch-Stoer algo- 
rithm for N-body simulations is that it is capable of keeping 
an upper limit on the local errors introduced due to taking 
finite time-steps by adaptively reducing the step size when 
interactions between the particles increase in strength. Cal- 
culations covered 1500 days, i.e. a total time span of avail- 
able transit observations. The results of simulations are pre- 
sented in Fig. [6] Our analysis allows us to exclude an addi- 
tional Earth-mass planet close to inner 1:3, 1:2, and 3:5 and 
outer 5:3, 2:1, and 3:1 mean-motion resonances (MMRs). 



5 LONG-TERM STABILITY OF THE SYSTEM 

In this section, we investigated the long-term gravitational 
influence of TrES-3b planet on a potential second planet in 
the system. Thus, we performed numerical simulation for 
studying the stability of orbits and checking their chaotic 




Figure 6. The upper-mass limit of a hypothetical additional 
planet that could perturb the orbital motion of TrES-3b as a 
function of ratio of orbital periods of transiting planet, P b , and 
the hypothetical perturber, P p . Orbits located in a grey area were 
found to be unstable due to close encounters of both planets. 



behavior using the method of maximum eccentricity (e.g. 
iDvorak et"aDl2003l ). 

We used long-term integration of small-mass (Earth- 
mass) planet orbits for inspecting the stability regions in 
the TrES-3b system. The time span of the integration was 
around 140000 revolutions of the planet around the star. 
For the integration of this system we again applied an ef- 
ficient variable-time-step algorithm (Bulirsch-Stoer integra- 
tion method) . The parameter e which controls the accuracy 
of the integration was set to 10 -8 , in this case. We have gen- 
erated 10 5 mass-less particles for representing small planets 
in this system. We assumed that semi-major axis ranges 
from 0.00625 au to 0.1 au. The inner border of the gener- 
ated system of small planets is around 1.5 times star radius. 
This value of inner border is approximately at the location 
of Roche limit for a Earth-like planet and includes also the 
region in the vicinity of 3:1 MMR (0.01 au for this case) in 
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our analysis. The step-size in semi-major axis was Aa = 
0.00025 au, in eccentricity Ae = 0.02 and A i = 5° in incli- 
nation. The upper limit of the grid in eccentricity was 0.5 
and 50° in inclination. We integrated the orbits of the small 
planets for the time-span of 500 years (about 140 000 rev- 
olutions of TrES-3b around the parent star). We obtained 
orbital evolution of each small planet in the system and also 
were able to study the stability regions at the end of the in- 
tegration using the method of maximum eccentricity, men- 
tioned above. The maximum eccentricity that the potential 



perturber of TrES-3b reached over the time of the integra- 
tion is plotted on the a — e (a — i) stability map in the Fig. [7] 
(Fig.®. 

Figure [7] shows the stability maps in a — e plane for se- 
lected values of inclinations: i = 5°, i = 20°, i = 35° and 
i = 50° (from top left to bottom right). Also the MMRs 
with TrES-3b planet are marked in these plots. We can see 
a stable region inside the 2:1 MMR for each value of in- 
clination. Outside the 2:1 MMR the gravitational influence 
of the TrES-3b planet is very strong and leads to deple- 
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tion of these regions. Only the MMRs 2:1, 3:2 and 1:1 are 
moderately populated but the population decreases with the 
increase of the inclination. For completeness, we note that 
5:3 (a « 0.016 au) and 3:5 (a « 0.032 au) MMRs are de- 
pleted and thus not stable for additional Earth-mass planet 
in their vicinity. In Fig.[S] we present the stability maps in 
a — i plane for several values of eccentricities: e = 0.001, 
e = 0.1, e = 0.3 and e = 0.5. One can see the stable region 
inside the 2:1 MMR with TrES-3b planet, which we found 
stable also in the a — e plane. Considering the region beyond 
the 1:3 MMR, the depletion of the planet population is not 
so strong than in the region between 2:1 and 1:3 MMR. This 
feature is in a good agreement with the weaker gravitational 
influence of the TrES-3b planet. 

Based on our results, we showed that the region inside 
the TrES-3b planet orbit (especially inside the 2:1 up to 3:1 
MMR) can be stable on longer timescales (hundreds of years 
or hundred thousands of revolutions of TrES-3b around the 
star) from the dynamical point of view. The region from 
0.015 au (near 2:1 MMR) to 0.05 au (near 1:3 MMR) was 
found unstable apart from moderately populated MMRs lo- 
cated in this area (see figures Fig. [7] and Fig.[8jl . The rela- 
tively small increase of population in the mentioned MMRs 
depends on the initial values of semi-major axis, eccentricity 
and inclination. The region beyond the 0.05 au was found to 
have a chaotic behavior and the depletion of the planet pop- 
ulation increases with increasing values of initial eccentricity 
as well as inclination. 



6 DISCUSSION AND CONCLUSION 

Based on the transit light curves obtained at several obser- 
vatories between May 2009 and September 2011 we rede- 
termined orbital parameters and the radius of the transit- 
ing planet TrES-3b. The best light curve (obtained at Calar 
Alto, Sep 2010) was used for light curve analysis, and the 
data from the other observatories were used for Tc deter- 
mination and TTV investigation. We used two independent 
solutions for parameters determination and finally, we con- 
cl uded that our values are consistent with previous r esults 
of ISozzetti et al.1 (|2009l ) and lChristiansen et all (120111 ). 

The aim of this present paper was also to discuss pos- 
sible presence of a second planet in the extrasolar system 
TrES-3. For this purpose we used our new 14 mid-transit 
times and t he individual determ inations from lSozzetti et"all 
(|2009h and IGibson et all i|2009l ). The resulting O -C dia- 
gram showed no significant deviation of data points from 
the linear ephemeris. In addition, we tried to search for a 
periodicity that could be caused by an additional body in 
the system. We can conclude that a strictly periodic TTV 
signal with the amplitude greater than 1 minute over a 4- 
year time span seems to b e unlikely. Th i s resu lt, together 
with refined assumptions of IGibson et al.l l|2009h allow us to 
put upper constraints on the mass of a potential perturb- 
ing planet. The additional Earth-mass planet close to inner 
3:1, 2:1, and 5:3 and outer 3:5, 1:2, and 1:3 MMRs can be 
excluded. 

Finally, we used the long-term integration of the the- 
oretical set of massless particles generated in TrES-3 sys- 
tem for studying the dynamical stability of potential sec- 
ond planet in the system (influenced by TrES-3b gravita- 



tion). From our analysis we found that the region inside the 
TrES-3b planet orbit (especially inside the 2:1 MMR) up 
to 3:1 MMR can be stable on longer timescales (hundreds 
of years or hundred thousands of revolutions of TrES-3b 
around the star) from the dynamical point of view. The 
region from 0.015 au (near 2:1 MMR) to 0.05 au (near 1:3 
MMR) was found to be unstable apart from moderately 
populated MMRs located in this area. The relatively small 
increase of population in these MMRs depends on the ini- 
tial values of semi-major axis, eccentricity and inclination. 
The region beyond 0.05 au was found to have a chaotic be- 
havior and depletion of the planet population increases with 
increasing values of initial eccentricity as well as inclination. 
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